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L ' invention cone erne, de fagon general e, les techniques 
de prospection par excitation sonique . 

Plus precisement , 1 1 invention concerne un precede de 
5 caracterisation geotechnique d'un fond sub- aqua ti que , tel 
qu'un fond sous-mar in, recouvert d'une nappe d'eau de 
hauteur totale de t ermine e , et comprenant une pluralite de 
strates formant entre elles des interfaces de separation, 
presentant chacune des parametres physiques propres et 

10 s 1 etendant a des profondeurs diverses en dessous d'une 
premiere strate qui forme une interface de separation 
avec l'eau, ce precede comprenant au mo ins une phase 
d 1 investigation, une phase de pre-traitement et une phase 
d ' analyse, la phase d ' investigation incluant elle-meme 

15 une operation d' excitation acoustique, mise en oeuvre en 
erne tt ant au moins vers le fond, depuis une source 
acoustique immergee, une onde acoustique incidente de 
signature f requentielle connue, une operation 

d * acquisition mise en ceuvre en produisant, au moyen d'une 

20 flute d'au moins quatre recepteurs immerges, des signaux 
de mesure respectifs resultant d'une detection d'ondes 
acoustiques respectives reflechies par le fond, et une 
operation de releve de donnees geometriques , la flute 
etant au moins approximativement alignee sur la source et 

25 distante du fond, les recepteurs etant espaces les uns 
des autres, et 1 ' onde incidente etant emise dans chaque 
fenetre d'une succession de fenetres temporelles 
d ' emission dis jointes , et presentant une variation 
temporelle de frequence a 1 « interieur de chaque fenetre 

30 d 1 emission. 

Un precede de ce type est par exemple decrit dans le 
brevet US 6 151 556. 
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25 



D-autres precedes sent egalement connus de 1'homme du 
metier, et notamment ceux que decrivent les documents de 
brevets EP - 0 215 703 et EP - 0 553 053. 

De fagon plus generale, les systemes acoustigues connus 
pour 1 ' investigation des fonds sont soit des systemes 
d'imagerie, soit des systemes sismigues. 



Les systemes d'imagerie classigues utilisent des ondes 
acoustigues de frequences elevees, typiquement 
superieures a 3 0kH 2 , emises vers le fond et permettant 
d'etablir une image de la surface du fond independamment 
du sous-sol en presence en se basant sur le mecanisme de 
15 • retro-diffusion. 

Les systemes sismigues usuels du type precedemment defini 
utilisent, en tant gu'emetteur, un boomer produisant un 
signal acoustique a tres basse frequence et presentant 
done l'avantage de penetrer profondement dans le sous- 



sol . 



Le signal issu des multiples reflexions et recueilli sur 
une flute d ' hydrophones et/ou de geophones tractee par un 
navire en synchronisms avec le boomer est ainsi 
representatif des couches prof ondes du sous-sol. 



du 



En revanche, les caracteristiques f requentielles _ 
signal emis a tres basse frequence sont mal controlees, 
30 de sorte que ces systemes ne fournissent que des 
indications qualitatives des strates profondes du sous- 
sol, a des profondeurs typiquement superieures a 10 
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metres, et avec des resolutions verticales typiquement 
superieures a 1 metre . 

Les sondeurs de sediment re cents presentent 1 ' avantage 
5 d'utiliser des signaux controles ("chirp" ou sources 
parametriques) qui insonifient le fond vertical ement a 
forte puissance. Toutefois, la reception des signaux, 
assuree par le meme transducteur , limite 

signif icativement les performances pour des raisons 
10 instrument ales (saturation electronique) et physiques 
(melange de phenomenes de retro-diffusion et de reflexion 
speculaire, ambigurte epaisseur/vitesse dans les 
couches) . 

15 Dans ce contexte, la presente invention vise a proposer 
un precede repondant au besoin, aujourd'hui non 
satisfait, de permettre la determination guantifiee des 
proprietes geotechniques du sous-sol d'un fond sub- 
aquatique sur les premiers metres, typiquement sur les 

20 cinq premiers metres au mo ins , et avec une resolution au 
plus egale a cinquante centimetres . 

A cette fin, le precede de l r invention, par ailleurs 
con forme a la definition generique qu'en donne le 

25 preambule ci-dessus, est essentiellement caracterise en 
ce que la phase d ■ investigation inclut une operation de 
displacement, c one omi tan t e avec les operations 

d' excitation et d ' acquisition, et mise en ceuvre en 
deplagant simultanement la source acoustique et la flute 

30 dans la nappe d'eau, en ce que la phase de pre- trait ement 
est mise en ceuvre en deduisant, a partir des signaux de 
mesure, de la signature de 1 ' onde acoustique incidente, 
et des donnees geometriques , des signaux de reflexion 
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ay ant atteint des recepteurs dif f erents a la suite d'une 
meme emission d r onde incidente, et en regroupant, sous 
forme de series de signaux de reflexion recales, des 
signaux de reflexion ayant atteint, sous des angles de 

5 rasance dif f erents, des recepteurs dif f erents en 
provenance d'une meme zone du fond a la suite de 
differentes emissions de l'onde incidente, et en ce que 
la phase d 1 analyse est mise en oeuvre en deduisant les 
parametres physiques de chaque strate du fond, a partir 

10 des series de signaux de reflexion recales, par 
exploitation de la dependance angulaire et f requentielle 
du coefficient de reflexion de I'onde incidente sur 
chaque interface, chaque signal de reflexion recale etant 
interprets comme une some d'arrivees retardees et 

15 attenuees de rayons propres correspondant a la reflexion 
de 1 ' onde incidente par les differentes strates de la 
zone du fond a laquelle correspond ce signal de reflexion 
recale ♦ 

20 

En pratique, chaque onde reflechie atteignant le centre 
de la flute dans une fenetre tempore lie de reception est 
avantageusement issue de la reflexion d'une onde 
incidente ayant atteint le fond sous un angle de rasance 
25 au plus egal a 70 degres et au moins egal a 10 degres , et 
de preference compris entre 2 5 et 65 degres. 

Chaque recepteur est par exemple separe de la source 
acoustique par une distance comprise entre environ 1 fois 
30 et environ 4 fois la hauteur d'eau moyenne sous les 
instruments, chaque recepteur pouvant par ailleurs etre 
espace d'un recepteur voisin par une distance comprise 
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entre environ 0.375 fois et environ 1 fois la hauteur 
moyenne suivant le nombre de ces recepteurs . 

Selon un autre mode d ' evaluation des distances, chaque 
5 recepteur est typiquement separe d'un recepteur voisin 
par une distance correspondant a une difference d ' angle 
de rasance au mo ins egale a 2 degres, et de preference au 
moins egale a 5 degres, des signaux de reflexion 
respectivement regus par ce recepteur et par le recepteur 
10 voisin. 

De preference, 1 1 onde incidente presente, dans chaque 
fenetre d 1 emission, au moins une composante de frequence 
relativement basse, dont la frequence est au moins egale 
15 a 10 0 Hz, cette onde incidente presentant utilement 
aussi, dans chaque fenetre - d l emission, au moins une 
composante de frequence relativement haute, dont la 
frequence est au plus egale a 3 kHz. 

20 Dans un mode de realisation possible de 1' invention, 
1 ' onde incidente peut ainsi presenter, dans chaque 
fenetre d 1 emission, une frequence continument variable 
entre une composante de frequence relativement basse et 
une composante de frequence relativement haute. 

25 

Pour obtenir a la fois une information suf f isamment riche 
et une separation satisf aisante des ondes reflechies, les 
fenetres temporelles d ' emission ont chacune par exemple 
une duree au moins egale a 0.1 seconde, et de preference 
30 au moins egale a 0.5 seconde, et sont repetees 
periodiquement avec une periode par exemple au moins 
egale a 0.5 seconde et de preference au moins egale a 1 
seconde . 
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Dans son mode de realisation prefere, l e precede de 
1' invention comprend, en amont de la phase d' analyse, une 
operation de determination de la profondeur d' immersion 
de la source acoustique et/ou de la profondeur 
d'immersion d'un point de la flute et/ou de la hauteur 



d ' eau . 



10 



15 



20 



25 



La phase d' analyse est par exemple raise en oeuvre en 
parcourant, de fagon iterative, une boucle incluant une 
operation de modelisation au cours de laquelle sont 
calcules, A partir de parametres physiques fictifs 
virtuellement attribues a la premiere strate au moins de 
la pluralite de strates, et pour les differents 
recepteurs, des signaux de reflexion recales fictifs 
respectifs issus de signaux de mesure virtuellement 
produits par ces recepteurs en presence de ces parametres 
fictifs, une operation de sous traction au cours de 
laquelle sont formes des ecarts constates entre les 
signaux de reflexion recales respectifs issus des 
signaux de mesure reellement produits par ces recepteurs 
et les signaux de reflexion recales fictifs 
correspondants, et une operation d - optimisation et de 
retroaction a 1- issue de laquelle les parametres 
Physiques fictifs sont revises en fonction des ecarts 
constates . 



30 



Le precede de !■ invention permet ainsi une 
caracterisation geotechnique de chaque strate d>une 
pluralite de strates d'un fond sub-aquatique au moyen de 
parametres physiques choisis dans 1 - ensemble de 
parametres qui comprend 1 < epaisseur de cette strate, la 
densite de cette strate, la vitesse du son dans cette 
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strate, 1 1 attenuation du son dans cette strate, la 
rugosite de cette strate, et un gradient de celerite du. 
son a l'interieur de cette strate. 

5 Dans son mode de realisation pref ere, le procede de 
1 1 invention permet d 1 identifier les differentes strates 
du sous-sol sub- aquatique avec des resolutions verticales 
de 1 ' ordre de dix centimetres . 

10 D 1 autres caracteristiques et avantages de 1 ' invention 
ressortiront clairement de la description qui en est 
faite ci-apres; a titre indicatif et nullement limitatif, 
en reference aux dessins annexes, dans lesquels : 

15 - la figure 1 est une vue schemata. que en perspective d'un 
equipement mettant en ceuvre le procede de 1 1 invention ; 

- la figure 2 est un diagramme schematique present ant de 
fagon approximativement chronologique les principales 

20 phases et operations du procede de 1 T invention ; 

- la figure 3, composee des schemas 3A a 3C, illustre de 
maniere intuitive et partielle 1 ' etape RECAL de recalage 
et de regroupement de signaux de reflexion mise en oeuvre 

25 pendant le deroulement du procede de 1' invention, cette 
etape conduisant notamment a regrouper des signaux de 
reflexion provenant d'une meme zone du fond; 

- la figure 4, composee des schemas 4A et 4B, complete 
30 1' illustration de 1' etape RECAL; 



- la figure 5, composee des schemas 5A et 5B, illustre 
parti el lement les mecanismes de propagation acoustique 



5 



20 



25 



30 
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multiple dans 1'eau, et leur exploitation dans 
1 ' invention; 



la figure 6 il lustre partielleiuent les mecanismes de 
propagation acoustigue multiple dans l e fond sub- 
aquatique; 




- les figures 7A a 7J sont des diagrammes representant, 
en fonction du temps t exprime en nombre d ' echantillons 

10 I- amplitude relative A de signaux de mesure 
respectivement produits par 10 receptees successifs de 
Plus en plus eloignes de la source acoustigue et separes 
les uns des autres par une distance importante, et 
mettant en evidence 1' influence des angles avec lesguels 

15 les ondes se ref lechissent au fond ; et 

- les figures 8A a 8 J sont des diagrammes representant, 
en fonction du temps t exprime en nombre d' echanti lions , 
1- amplitude relative A de signaux de mesure 
respectivement produits par 5 receptees successifs de 
Plus en plus eloignes de la source acoustigue et pour une 
modulation de basses frequences a gauche et pour une 
modulation de frequences hautes a droite, mettant en 
evidence 1 < influence des frequences utilisees. 

Notations 

Le lecteur trouvera ci-apres guelques-uns des symboles 
utilises pour identifier les differentes entites 
impliquees dans la mise en oeuvre de 1' invention : 

"i- est un indice compteur utilise pour numeroter chague 
hydrophone en le distinguant ainsi des autres 
hydrophones ; 
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"j n est un indice compteur utilise pour numeroter chaque 
strate du fond sub-aquatique en la distinguant ainsi des 
autres strates du fond; 

5 

M J" est un nombre entier designant le nombre total de 
strates observables ; 

"k" est un indice compteur utilise pour numeroter chaque 
10 rayon propre d'une onde acoustique en le distinguant des 
autres rayons propres ; 

"1" est un indice compteur utilise pour numeroter chaque 
emission d'une onde acoustique incidente (encore appelee 
15 "tir") en la distinguant ainsi des autres emissions d'une 
onde acoustique; 

"n" est un indice compteur utilise pour numeroter chaque 
boucle d 1 iteration de la phase d' analyse en la 
20 distinguant ainsi des autres boucles d' iteration; 

H e" est un indice attribue a une valeur finale exacte 
d'un parametre; 

s(t) designe 1 1 onde acoustique incidente; 

25 

si(t,l) designe le signal de me sure regu sur le i® me 
hydrophone a la suite du l ierne tir; 

ei ( t , 1 ) designe 1 ' enveloppe de 1 1 intercorrelation entre 
30 1 1 onde acoustique incidente s ( t } et le signal de mesure 
sift,!) ; 



15 
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rfi(t,l) designe le signal de reflexion associe au signal 
de mesure Si (t,l) par application a l'enveloppe ei (t,l) 
d'une fenetre de reception; 

5 Rfi(t,B) designe un signal de reflexion recale, issu du 
signal de reflexion rf i( t,l) par une etape de recalage 
RECAL et associe au fond B; 

RfSi<t, B) designe un signal de reflexion recale fictif, 
10 obtenu par simulation; 

Rf Si k (t, B ) designe le k"- ray on propre du signal de 
reflexion recale Rfi(t,B); et 



0i k (j) designe 1' angle de rasance sur la j i&ne strate du 
k rayon propre du signal de reflexion recueilli sur le 



1 



hydrophone . 



Comme indigue precedemment , 1 -invention concerne (figure 
20 1) un precede de caracterisation geotechnigue d'un fond 
sub-aquatigue B, tel qu'un fond sous-xaarin, reconvert 
d'une nappe d ■ eau de hauteur totale H determinee. 

Ce fond B comprend une pluralite de strates telles que 
25 Fl, F2, etc. (figure 6), generiquement notees Pj et 
presentant chacune des parametres physiques propres. 

La premiere strate PI forme une interface avec 1'eau, les 
autres strates s « etendant en dessous de cette premiere 
strate a des profondeurs diverses . 



30 



Le precede de 1 • invention vise, pour chacune des strates 
Fj d'un ensemble de strates comprenant typiquement, mais 
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non limitativement , entre deux et cinq strates, a 
determiner plusieurs parametres distinctif s de cette 
strate , tels que l'epaisseur hj de cette strate, sa 
densite p j , la vitesse du son Cj dans cette strate, et 
5 plus precisement la vitesse C P j pour des ondes de 
compression et la vitesse C S j pour des ondes de 
cisaillement , 1 ' attenuation du son ofj dans cette strate, 
et plus precisement 1 x attenuation a P j pour des ondes de 
compression et 1 ' attenuation a sj pour des ondes de 
10 cisaillement, ainsi, eventuellement , que la rugosite dans 
cette strate, ou un gradient de celerite du son. 

Le precede de 1 1 invention comprend essential lement une 
phase d' investigation INVESTIG, une phase de pre- 
15 traitements PRE_TRAIT et une phase d 1 analyse AJNTA. 

La phase d 1 investigation IWESTIG inclut elle-meme une 
operation de deplacement DEPL, une operation d 1 excitation 
acoustique EXCIT, une operation d 1 acquisition ACQUI, et 
20 une operation RE LEV de r el eve de donnees non acoustiques, 
ces operations etant raises en ceuvre de fagon global ement 
simultanee. 

L ' operation de deplacement DEPL consiste a deplacer dans 
la nappe d'eau, au moyen d'un navire N et a vitesse 
25 reduite ideal ement constante, une source acoustique 1 et 
une flute 2 au-dessus du fond B, la flute 2 comprenant 
par exemple trois, quatre, ou cinq recepteurs acoustiques 
21 a 24, encore appeles hydrophones. 

30 L l operation RELEV de releve de donnees non acoustiques 
comprend elle-meme une operation CAPT de capture de 
donnees environnementales et une operation de reperage 
REPER . 
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L ' operation CAPT de capture de donnees environnementales 
consiste a enregistrer a chaque instant t les profondeurs 
de la source acoustique 1 et d'un point de la flute 2, 
5 ainsi que la hauteur d'eau H . 

L 1 operation de reperage RE PER consiste a enregistrer, a 
chaque instant t , les coordonnees X et Y du point du 
navire N. 

10 

L 1 operation EXCIT d' excitation acoustique est mise en 
oeuvre en emettant de fagon aussi isotrope que possible, 
notamment vers le fond B et depuis la source acoustique 
immergee 1 , 1 1 onde acoustique incidente s (t ) f de 
15 signature f requentielle connue. 

L 1 operation d ' acquisition ACQUI est, quant a elle, mise 
en osuvre en produisant, au moyen des recepteurs immerges 
21 a 24, des signaux de mesure respectifs si ( t , 1) , 
20 s2{t,l), s3 (t, 1) et s4(t/l) ( generiguement notes Si ( t } ) 
en reponse a une detection d'ondes acoustiques 
respectives propagee dans 1'eau, et en enregistrant ces 
signaux en fonction du temps t. 

25 Comme le montre la figure 1, la source 1 et la flute 2 
restent a distance du fond B et se trouvent 
respectivement immergees a des profondeurs 2s et Zr, 
eventual lement egales, la flute 2 etant au moins 
approximativement alignee sur la source 1 et au moins 

30 approximativement parallele au fond B, 

Les profondeurs d 1 immersion de la source 1 et des 
recepteurs 21 a 24 sent par example enregis trees par des 
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capteurs de press ion, la profondeur d'eau H etant par 
exemple enregistree a l'aide d 1 un sondeur de bathymetrie. 

La hauteur totale de la nappe d'eau au-dessus du fond B 
5 etant notee H, la source 1 et la flute 2 se trouvent 
ainsi globalement , par rapport au fond B, a une hauteur 
moyenne Hm typiquement comprise entre 10 metres et 2 0 
metres . 

10 Les recepteurs 21 a 24 sont espaces les uns des autres, 
chaque recepteur etant separe d'un recepteur voisin par 
une distance d qui peut etre constante ou evolutive le 
long de la flute 2 . 

15 La flute 2 est consideree comme delimitee entre les deux 
recepteurs extremes 21 et 24, qui sont done equidistants 
du centre 20 de cette flute. On note Di, la distance 
horizontals separant 1'emetteur du i 6rne hydrophone et Zri 
la profondeur du i^ me hydrophone. 

20 

La phase d l analyse ANA, qui est destinee a permettre 
1 1 identification des parametres physiques tels que h j , 
pj , Cj, (Xj de chaque strate Fj du fond B , exploite a cette 
fin la signature f requentielle de 1 ' onde acoustique 
25 incidente s{t), et les signaux de mesure Si(t,l) des 
dif f erents recepteurs 21 a 24. 

Le precede de 1 1 invention concentre d'abord son analyse 
sur les signaux Si(t) propages entre la source 1 et les 
30 recepteurs 21 a 24. 



Pour ce faire, 1 ' onde incidente s ( t) , aussi isotrope que 
possible t est emise par la source 1, de directivite 
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connue, dans des fenetres temporelles d' emission 
successives dis jointes, et le signal de mesure Si(t,l) de 
chaque recepteur est lui-meme analyse dans des fenetres 
temporelles de reception successives dis jointes , au cours 
5 desquelles sont detectes les signaux propages . 

Dans ces conditions, chaque signal regu Si(t) est 
constitue de repliques attenuees et retardees du signal 
emis s ( t ) pour le tir note Tl , la structure du signal 
10 regu Si(t) etant ainsi fonction de la geometrie du trajet 
multiple que 1 • onde acoustique parcourt dans 1 1 eau 
(figure 5) . 

Ce trajet multiple comprend principalement un trajet 
15 direct de 1 < onde acoustique depuis la source 1 jusqu 1 a 
chacun des recepteurs 21 a 24, un trajet dit "reflechi 
fond" de 1 1 onde acoustique depuis la source 1 jusqu * a 
chacun des recepteurs 21 a 24 via reflexion sur le fond, 
1 1 onde acoustique ayant parcouru ce trajet dormant 
20 naissance au signal de reflexion rfi(t,l), et un trajet 
dit "reflechi surface" de 1 ' onde acoustique depuis la 
source 1 jusqu 1 a chacun des recepteurs 21 a 24 apres 
reflexion a la surface de l l eau. 

25 Dans le cas particulier des fonds a faible hauteur d f eau 
H ou "petits fonds", on utilisera avantageu semen t le 
trajet de reflechi surface, voire le trajet dit de 
"reflechi fond- surface 1 ' impliquant une reflexion de 
1 • onde acoustique a la fois sur le fond et sur la surface 

30 entre son emission par la source 1 et sa detection par 
1 ' un des recepteurs 21 a 24. 
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La mesure des temps d'arrivee respectifs de ces trajets 
multiples dans 1 1 eau permet d'une part de determiner avec 
precision la geometrie du dispositif a 1 ' instant de la 
mesure et d' autre part d'etablir des fenetres de 
5 reception a 1 1 interieur desquelles se concentrera 
1 1 analyse des ref lechis fond. 

Plus precisement , le procede de 1 1 invention concentre son 
analyse sur les signaux de mesure issus d'ondes 

10 ref lechies une seule fois par le meme fond B , Comme 
indique sur la figure 3, dont les differents schemas 3A a 
3C correspondent a des tirs differents, un ordonnancement 
de plusieurs tirs consecutifs permet alors de recueillir 
les signaux notes Rf i (t,B) issus de la reflexion de 

15 1 1 onde acoustique incidente sur une meme zone de fond B 
et sous autant d 1 angles de rasance differents qu ' il 
existe de recepteurs tels que 21 a 24. 

Comme le montre la figure 4, les differents signaux de 
20 reflexion rfi{t,l) recueillis, lors d'un tir 1, sur 
chacun des recepteurs 21 a 24 dans une fen§tre 
d ' observation correspondant a la reception d'un echo 
unique sur le fond doivent done §tre respect ivement 
associes a differents signaux de reflexion r£i(t,l'), 
25 rf i (t, 1 1 1 ) , rf ± ( t , 1 ' ? ' ) etc., recueillis, lors d'autres 
tirs 1', 1 ' ' , I 1 1 ' , etc. sur chacun des recepteurs 21 a 
24 pour former les signaux dits de "reflechi fond" ou de 
reflexion recales Rfi(t,B), cette operation etant 
effectuee a une etape de recalage REGAL qui sera 
30 detaillee ulterieurement. 

En realite, chaque signal de reflechi fond Rf i (t / B) est 
lui-meme composite en ce sens qu'il est forme (figure 6) 
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par la propagation, suivant un trajet multiple, de 1 1 on.de 
acoustique incidente dans le fond B / dormant lieu a des 
reflexions de cette onde incidente sur les differentes 
strates Fj du sous-sol sufo-aquatique . 

5 

Dans ces conditions, chaque signal de ref lechi fond 
Rf i ( t , B) peut etre decrit sous la forme d'un ensemble de 
rayons propres notes Rf Si k (t , B) , 1 ' indice k se rapportant 
au k ieme rayon propre associe a 1 1 onde incidente Rfi(t,B), 

10 Le rapport d' amplitude entre le signal temporel incident 
et ref lechi a 1 7 interface Fj est directement 
proportionnel a la transformee de Fourier du coefficient 
de reflexion Rj (f ,6i k { j } ) , 1' indice k se rapportant ici a 
1 ' angle 9 t k ( j ) avec 1 equal le k xeme rayon propre vient se 

15 ref lechir sur 1' interface F j . 

Or, ce coefficient de reflexion Rj { f , 0i k ( j ) ) , dans le 
cadre d'une modelisation de Rayleigh, est relie aux 
parametres geotechniques des differentes strates Fj-1 Fj 
20 du sous-sol. 

Les signaux ainsi produits et exploites portent la 
signature de 1 ' environnement rencontre par les ondes 
detectees , et contiennent done une information 
25 suf f isamment riche et pertinente du fond B pour permettre 
d'en extraire sans ambiguite les parametres geotechniques 
des differentes strates Fj du sous-sol. 

L ' echo d'une onde incidente parvenant au fond B sous un 
30 angle de rasance 6 tres rasant n ' est pratiquement 
sensible qu 1 a 1 ' attenuation acoustique de la strate Fl . 
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L ' echo d ' une onde incidente parvenant au fond B sous une 
incidence perpendiculaire au fond B n ' est prat iquement 
sensible qu'aux changemen t s d' impedance acoustique et 
presente notairiment le defaut de ne pas pouvoir permettre 
5 sans ambigulte la determination du couple epaisseur- 
vitesse de chacune des strates . 

De f agon generale , il est done preferable de recuperer 
des signaux ayant subi une reflexion a des angles de 
10 rasance assez faibles sur chacune des strates Fj 
compos ant le fond. 

Ainsi, des angles de rasance moderes sent neeessaires 
pour decrire la premiere strate Fl et des angles de 

15 rasance plus import ants sont neeessaires pour decrire les 
couches plus prof onde s , F2 , etc. II est done apparu tres 
opportun de faire en sorte, pour 1' invent ion, que chague 
onde ref lechie atteignant l'un des recepteurs 21 a 24 
dans une fenetre temporelle de reception soit issue de la 

20 reflexion d'une onde incidente ayant atteint le fond B 
sous un angle de rasance 8k au plus egal a 70 degres et 
au moins egal a 10 degres. 

Bien que les figures 5 et 6 soient a la fois totalement 
25 schematiques et incompletes compte tenu de leur mode de 
presentation simplif ie et purement spatial, elles 
illustrent, de fagon intuitive, ce que sont les trajets 
multiples de 1 1 onde incidente, respect ivement dans 1 1 eau 
et dans le fond. 

30 

L ' importance ci-dessus discutee des angles d 7 emission vis 
a vis du ref lechi fond peut etre observee sur les courbes 
des figures 7A a 70* , correspondant respectivement aux 
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ondes recueillies sur des recepteurs disposes a des 
distances correspondant a des angles de rasance au fond 
decroissant de 72 a 25 degres . 

Ainsi, il apparait que la courbe 7B, qui correspond a des 
angles de rasance forts, apporte peu d' information 
relativement a la courbe 7 A puisque pour ces angles, le 
coefficient de reflexion de chaque strate est 
pratiquement independant de 1' angle de rasance. De meme, 
il apparait que les courbes 71 et 7J qui correspondent a 
des angles de rasance faibles n'apportent plus 
d' information puisque les ondes ne penetrent pas dans le 
fond, inversement, les figures 7B a 71 qui correspondent 
a un balayage angulaire pertinent, presentent une forte 
variability permettant de retrouver sans ambiguite les 
caracteristiques du fond. 

De meme, alors que la signature d'une onde incidente tres 
basse frequence est difficile a maitriser, et qu ■ une onde 
a haute frequence presente une capacite de penetration 
nulla dans les strates prof ondes, il est apparu tres 
opportun, pour !• invention, d'utiliser une onde incidente 
presentant une composante de frequence relativement basse 
au moins egale a 100 H 2 , et une composante de frequence 
relativement haute au plus egale a 8 kHz, cette onde 
incidente presentant par example, dans chaque fenetre 
d- emission, une frequence contingent variable entre 100 
Hz et 8 kHz. 

La relation liant la profondeur de penetration d'une onde 
acoustique incidente a la frequence de cette onde est 
clairement illustree par les figures 8A a 8J, qui 
representent, pour 5 angles de rasance 9i differents de 
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1 1 onde incidente compris entre 25 et 70 degres , la 
reponse Rfi(t,B) du signal temp or el reflechi au fond, 
pour un signal d' emission compris entre 400 et 3200 Hz 
(figures 8A a 8E) et compris entre 5400 et 8200 Hz 
5 (figure 8F a 8 J) . Bien que presentant la rueme resolution 
temporelle que les signaux basse frequence, les signaux 
haute frequence sont sensiblement appauvris par 
1 ' attenuation des ondes dans les diverses strates. 

10 Ainsi , grace a un choix approprie de 1 ' angle de rasance 
et grace a une variation de la frequence acoustigue 
d* excitation, les parametres physiques hj , pj , C j , cxj de 
chaque strate Fj peuvent etre deduits de 1 * onde incidente 
s(t) et des signaux de mesure Si(t,l) par 1 1 intermediaire 

15 d'un model e physique faisant intervenir les coefficients 
de reflexion acoustigue Rj (6i k ( j ) ) des strates Fj du fond 
B vu par les rayons propres Rfs i k (t,B) / leur attenuation 
aj, ainsi que leur dependance f requentielle . 

20 Autrement dit, le modele physique employe permet de tirer 
parti des dependances angulaires du coefficient de 
reflexion et de la dependance frequentielle de 
1 ' attenuation dans les strates. 

25 En pratique, il est souvent utile de faire en sorte que 
chaque onde ref lechie atteignant le centre 20 de la flute 
2 dans une f enetre temporelle de reception soit issue de 
la reflexion d'une onde incidente ayant atteint le fond 
sous un angle de rasance 9j. au plus egal a 70 degres et 

30 au moins egal a 10 degres, 1 ' angle de rasance Qi pouvant 
meme le plus souvent etre idealement choisi dans une 
plage angulaire comprise entre 25 et 65 degres. 
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Chaque recepteur 21 a 24 est alors separe de la source 
acoustique 1 par une distance D comprise entre environ 1 
fois et environ 4 fois la hauteur moyenne Hm, la valeur D 
= Hm concernant essentiellement le recepteur 21 le plus 
5 proche de la source 1, et la valeur D = 4 . Hm concernant 
plut6t le recepteur 24 le plus eloigne de la source 1. 

Dans ces conditions, chaque recepteur, tel que 21 a 23, 
est espace d'un recepteur voisin, tel que 22 a 24, par 
10 une distance d permettant une discretisation la plus 
reguliere possible de 1' angle de rasance. 

En pratique, les fenetres temporelles d' emission ont 
chacune une duree au moins egale a 0.1 seconde, et de 
preference plutot au moins egale a 0 . 5 secondes . De f agon 
generale, la duree choisie est un compromis entre un 
signal court minimisant les effets du mouvement du navire 
et un signal plus long augmentant le rapport signal a 
bruit . 

Ces fenetres d« emission sont repetees periodiquement avec 
une periode dependant de leur duree individuelle , au 
moins egale a 0.5 secondes, et de preference plutot au 
moins egale a 1 seconde. De facon generale, la duree du 
silence entre deux tirs est un compromis entre une duree 
elevee permettant une separation satisf aisante des 
signaux recus et une duree faible permettant une plus 
grande resolution au sol. 

Dans le cas prefere, les signaux de mesure sont analyses 
par tranche de 2 secondes et traites par tir, un tir 
correspondant a un jeu de signaux forme d'un signal emis 
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et des signaux de mesure correspondants regus par les 
recepteurs respectifs 21 a 24. 

Entre la phase d ' investigation INVESTIG telle que 
precedemruent decrite, et la phase d* analyse ANA qui sera 
davantage detaillee ulteri eurement , le procede de 
1 ' invention comprend une phase de pre-traitement notee 
PRE_TRAIT . 

Cette phase de pre-traitement PRE_TRAIT comprend une 
operation de filtrage FIL_COR, une operation de calcul 
geometrique INVER_JGE0M, et une etape de recalage REGAL. 
Au cours de 1 ' operation FIL^COR, chaque signal de mesure 
si(t) est filtre pour normal iser ce signal en tenant 
compte de la signature frequent ielle de 1 ' onde acoustique 
incident e s (t) . 

Au cours de 1 1 operation 1WER_GE0M, la geometrie de 
1 ' instrumentation lors du tir est recalculee et les 
fenetres temporelles de reception sont identif iees pour 
pouvoir extraire les signaux de reflexion rfi(t,l)« 

Enf in, au cours de 1 ' etape REGAL , chaque signal de mesure 
extrait est associe a une position geographique du fond 
reperee par ses coordonnees X et Y , et recale de f agon a 
tenir compte du deplacement du navire entre chaque tir et 
le tir suivant, 1 » etape REGAL produisant ainsi un jeu de 
signaux de reflexion recales Rfi(t,B) ayant vu une meme 
zone de fond B lors de tirs differents et sous des angles 
de rasance differents. 

Le filtrage effectue lors de 1 1 operation FIL_COR, dit 
"adapte" , est obtenu par estimation de 1 ' enveloppe de 
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l'intercorrelation entre le signal emis represents par 
l'onde acoustique incidente s(t), et chaque signal de 
mesure Si (t , 1) . 

L ' operation de filtrage consiste d'abord a effectuer la 
5 correlation : 
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c,. (/, 0 = js(r).s t (t - x, l)dr , 

puis a effectuer une extraction d'enveloppe du signal 
Ci(t,l) qui se realise par exemple a 1 ' aide d' une 'phase 
de demodulation selon C< (/, , ou f c est la frequence 
centrale du signal emis et suivie d'un filtrage passe-bas 
a la frequence de coupure F coup > (Fmax-Fmin) /2 pour 
obtenir le signal ei (t, 1) . 

Le calcul de la geometrie instrumentale exacte effectue 
lors de 1' operation INVER_GEOM est base sur 
1' identification des ondes reflechies arrivant 
successivement sur les differents recepteurs 21 a 24. 

Cette identification est obtenue par exemple a partir 
d'une detection de maximum a 1 ' interieur de fenetres 
temporelles qui sont definies a partir de la geometrie a 
priori de 1 7 instrumentation . 

On compare les temps mesures (^) m correspondant a cette 
detection aux temps (xe^) . issus d'un modele de rayons 
propres en milieu stratifie fournissant 1' expression de 
la reponse impulsionnelle theorique : 

RreX0 = ±4*.S(t-re*) Equation 1 
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ou N represente le nombre de trajets acoustiques pouvant 
relier un recepteur a 1 ' emetteur , ou t represente la 
variable temps, ou xe^ represente le retard que 
presente, par rapport a 1 1 emission de 1 ' onde incidente, 
la reception d'une onde ayant suivi le k ieme trajet, ou 5 
represente la fonction de Dirac, et ou A ± k represente 
1' amplitude relative, rapportee a celle de 1 1 onde 
incidente emise, de 1 • onde ayant suivi le k idme trajet et 
parvenant au recepteur avec un retard xei k , la notation Te 
rappelant qu'il s'agit des retards associes au trajet 
multiple de 1 1 onde acoustique dans l'eau. 

Les retards associes au trajet direct et au reflechi fond 
sont successivement donnes par : 



*e, =Jti et Te i = 7, 

^moy ^mqy 

ou Cmoy est la vitesse moyenne du son dans la colonne 
d'eau traversee par 1 1 onde ayant suivi ce trajet et dans 
le cas particulier des grands fonds ou les reflexions en 
surface ne sont pas exploitees . 

Les instants theoriques prevus pour les differents rayons 
propres calcules sont compares aux instants d 1 arrivees 
des ondes reellement detectees par les differents 
recepteurs 21 a 24 . 

En cas de difference non negligeable, une recherche de la 
geometrie exacte peut etre effectuee par example par un 
algorithme de recuit simule a 2.N variables (Di et Sri) 
en minimisant une fonction de cout exprimee a partir des 
differences entre les instants theoriques associes aux 
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rayons propres et les instants de detection des ondes 
recueillies sur les N differents recepteurs 21 a 24. 

Dans le cas particulier des petits fonds, d'autres 
arrivees telles que le reflechi -surface et meme les 
reflechis surface-fond ou fond-surface sent detectes et 
permettent de s'affranchir des mesures prealables de H et 
Zs, grandeurs qui seront directement recherchees par 
l'algorithme d' inversion geometrique decrit ci-dessus. 

A partir de 1 ' information alors disponible, qui- est 
constitute, pour chaque tir et chaque recepteur par 
le signal Si (t) issu de ce recepteur et par la geometrie 
du systeme, il convient de selectionner la plage 
temporelle correspondant a la fenetre de reception dans 
laquelle le signal regu sera analyse, e'est-a-dire la 
plage temporelle de reception de la seule contribution du 
signal de reflexion rfi(t,l). 

Cette selection est par exemple effectuee, dans le cas 
particulier des petits fonds, par un programme qui, pour 
chaque enveloppe du signal ei (t,l) issu du i^ e recepteur 
selectionne la fenetre temporelle definie par : 

>f, (?, /) = e, {ref - Nt e : re* - Nt, , I) 

ou est un separateur symbolique place entre la borne 

temporelle inferieure et la borne temporelle superieure 
de la fen§tre de reception selectionnee , ou Tei 3 
represente le retard associe au reflechi fond, ou x ei 4 
represente le retard associe au premier reflechi surface- 
fond, ou N est un parametre dependant de la nature du 




1 er depot 

25 



signal emis, et ou t e est la periode d' echantillonnage de 
chague recepteur, les indices 3 et 4 etant lies a un cas 
particulier des petits fonds dans leguel 1 ' ordre 
d'arrivee des rayons propres serait successivement le 
trajet direct, le ref lechi-surf ace, le ref lechi-f ond et 
le ref lechi surface-fond. 
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L' operation de recalage RECAL consiste a utiliser les 
donnees de position (X,Y) avec les positions des 
instruments associe au tir correspondant de facon a 
calculer les positions B^^Y) du fond vu par 
1' hydrophone i lors du tir 1. En repetant 1' operation 
pour plusieurs tirs consecutifs, il devient possible de 
reconstituer le jeu de donnees Rfi(t,B) constitue des i 
reflechis fonds ayant vu le fond B(X,Y) sous les 
differents angles de rasance 9i . 
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A titre d'exemple, pour une configuration definie par une 
distance de 10m entre la source 1 et le premier 
hydrophone 21 de la flute, pour une flftte a 4 hydrophones 
separes les uns des autres de . 10m, horizontale et situee 
a 10m au dessus du fond, pour un navire se deplacant a la 
Vitesse reguliere de 5m/s et avec une duree de 2 secondes 
entre tirs successifs, on appliguera successivement un 
decalage de 1 tir pour le 2 i&ne hydrophone 22, de 2 tirs 
pour le 3 1&ne hydrophone 23, et de 3 tirs pour le 4 i ^ me 
hydrophone 24 comme rappele dans le Tableau 1 suivant . 



KT ro de tir 


Signaux 
ext raits 


Decal . 


Recalage 


Tir 1 


rfi(t,l) 


0 


rf^t.l) -> Rf^t.B) 


Tir 1+1 


rfift.l+l) 


1 


rf 2 (t, 1+1) ~> Rf 2 (t,B) 


Tir 1+2 


rfi(t,l+2) 


2 


rf 3 (t / l+2)-> Rf 3 (t,B) 
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Tir 1+3 


rfi(t,l+3) 


3 


rf 4 (t,l+3)-> R£ 4 (t,B) 




Tableau 1: Principe de l'operation de recalage des tirs 

Une fois identifiee et extraite, pour chaque recepteur, 
l'onde detectee ayant subi une reflexion unique et par le 
fond B, la phase d' analyse ANA est raise en ceuvre pour 
permettre de deduire de cette onde les parametres 
geotechniques du fond au moyen d'un modele fonde sur 
1' expression du coefficient de reflexion de Rayleigh. 



Comme evoque precedemment , la modelisation des 
interactions des ondes acoustiques avec le fond et ses 
differentes strates est basee sur la representation de 
chaque strate sur une base de plusieurs parametres 
15 incluant l-epaisseur h de cette strate, la densite p {par 
exemple en g/cm 3 ) ; la vitesse du son C p/S (par exemple en 
m/s) et 1 ' attenuation du son a p/s (par exemple en dB/A) ; 
les indices p et s se rapportant respectivement aux ondes 
de compression et aux ondes de cisaillement . 



Dans les couches sedimentaires , la rugosite ou encore un 
gradient de celerite peuvent dans certains cas etre 
ajoutes a ces principaux parametres physiques. 



Partant du fait que le reflechi fond se decompose lui- 
meme comme une some de rayons propres reflechis par les 
differentes strates du fond, on peut egalement en ecrire 
sa reponse impulsionnelle selon : 
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Tbi k se rapportant au fait qu'il s'agit du trajet multiple 
de l'onde acoustique dans le fond. 

Considerant maintenant E(t) , 1 ' autocorrelation de 
l'enveloppe du signal emis, on obtient Rfi(t) par : 



Rfs ,(*, B) = E{t) ® £ 4.S(t - < ) = f; 4-E(t-Tb*) = f]rfs i k (t,B) 
k=i M w 

l'indice s precisant qu'il s'agit du reflechi fond 
simule. 



Concernant la determination de ces rayons propres, il 
s'agit en particulier de calculer les k angles de tirs 
(angle de rasance au niveau de la source 1) des k rayons 
propres recherches . 

On part de 1 ' angle Q ± ° defini par : 

On effectue alors une recherche dichotomique de 1 ' angle 
de tir 9i k (0) entre 9/ et O^+flUx en utilisant les 
relations : 



c °sg/C/> cos^C/4-1) 

~~ entre les strates Fj et Fj+1, en 

notant 0 d k (j) 1 'angle dans la couche Fj du rayon partant 
avec un angle de tir 0 i k (O) / et 
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D __ (2.H-Zr i -Zs) y 2.N(J,k).h(j) 
tan<?*(0) n tane, k (f) 

en considerant ici le k lfeme rayon propre &u reflechi fond 
associe a la geometrie du capteur i comme le rayon propre 
5 se ref lechissant au plus profond sur la j ldme strate, et 
en notant N(j,k) le nombre de reflexions du rayon k dans 
la couche F j . 



Un test doit prealablement etre effectue pour s 'assurer 
10 de 1 ' existence de ce rayon N (1, k) . 

Dans le cas ou Ai est une fonction de la frequence, on 
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Af(t) = FFT-\4(f)) 

Dans le cas particulier d' interfaces planes et sans 
reflexion multiple dans aucune couche, 



A k tn = (.-j^r~ )* ~ -e ^ .vj (/, 0, k (J - 1)). fl (i - rf (/> 0 " U - i))) 

ou dj est la distance parcourue dans la strate F j , 
tj=dj/Cj le temps correspondant, QCj 1 ' attenuation du son 
25 dans la couche exprimee en Neper /m, et Rj ( f , 9i k ( j ) ) le 
coefficient de reflexion a la j leme interface. 



Le coefficient de reflexion de Rayleigh associe a une 
interface f luide-f luide est defini par : 
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Pj+\- c j+\ P}* c j 

. sm(g,*(7 + l)) sin(g f A g)) 

sin(^0 + l)) sin(^(y)) 
ou 1 1 on utilise les celerites complexes : 



1 + - 



2.^.8,686 



Ce xnodele etant pose en amont de toute investigation, la 
phase d ! analyse ANA peut l'utiliser pour extraire les 
parametres geo techniques des differentes strates du sous- 
10 sol sub-aquatique a partir des signaux de reflexion 
recales et corriges Rf i (t / B) ofotenus a la fin de la phase 
PRE TRAIT. 



Pour ce faire, la phase d 1 analyse ANA utilise par exemple 
15 une technique d 1 appariement de formes connue de 1 ' homme 
de metier, auquel cas elle consiste a parcourir, de f agon 
iterative, une boucle incluant une operation de 
modelisation MODEL, une operation de sous traction SOUST, 
et une operation d ' optimisation et de retroaction 
20 OPTI_RETR . 



On rappelle que, par convention, l'indice "n" sera 
utilise ci-apres pour designer de fagon generique le 
n ombre courant de fois que la boucle d ' iteration a ete 
25 parcourue depuis le debut de la phase ANA, et pour 
reperer toute grandeur obtenue au n^ me pas du parcours 
iteratif de cette boucle. 
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La model isation MODEL, coxisiste a calculer, a chaque pas 
n, les signaux de reflexion recales fictifs respectifs 
Rf i n (t, B n ) que produiraient les recepteurs 21 a 24 si le 
fond B etait constitue par un fond note B n forme de 
5 strates Fj definies par des parametres physiques fictifs 
h n j ( p n j , C n j , ot n j virtuellement attribues a ces strates. 

Les parametres physiques fictifs initiaux h oj , p oj , C 0 j , 
a oj sont par exemple choisis sur la base d 1 indications 

10 geologiques generales sur le site a explorer, ou a defaut 
au hasard dans des f ourchettes de valeurs plausibles ou 
simp lament possibles, les valeurs prises par les 
parametres h nj/ p nj , C n j , a nj pour les valeurs non nulles 
de n etant definies de proche en proche par le parcours 

15 repete de la boucle d' iteration. 

A ces parametres peut etre ajoute le nombre J n defini 
comme le nombre total de strates supposees , lors de 
1 ' execution du n ldme pas de la boucle d 1 iteration, avoir 
20 ete detectees, ce nombre evoluant jusqu'a 1 ' obtention, a 
1' issue de 1 1 iteration, du nombre total exact J e de 
strates reellement detectees . 

Pour 1' obtention des parametres exacts, la phase 
25 d 1 analyse ANA comprend tout d'abord une operation de 
soustraction SOUST au cours de laquelle sont calcules les 
ecarts A ni que presentent, par rapport aux signaux de 
reflexion recales fictifs Rf si n ( t , B n ) prealablement 
calcules, les signaux de reflexion recales f litres et 
30 corriges Rfi(t,B e ). 

Puis, la phase d ' analyse ANA comprend une operation 
d ' optimisation et de retroaction OPTI_JRETR a 1' issue de 
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laquelle les parametres physiques fictifs sont revises en 
f onction des ecarts constates . 

Pour ce f aire, le procede de 1 1 invention calcule, a 
5 chaque pas d ' iteration n, une fonction de cout donnee par 
exemple par: 



10 avec A ni = RfSi n {t, B n ) - Rfi(t,B e ) • 

Cette fonction de cout est par exemple minimi see selon un 
algorithms de recuit simule pour produire une reponse 
impulsionnelle affinee Rfs i n+1 (t / B n+ i) definissant un 
15 nouveau fond B n+i utilise dans un nouveau pas d* iteration 
de la phase ANA, le processus convergeant vers la 
production de parametres exacts h e j/ p e j/ C e j / cx ej des 
dif ferentes strates j . 

20 Le tableau ci-dessous donne un exemple de cas concret 
traite . 



Geometrie : 

H=25m, Zs=7m / Zr i =8m / D i =12m / di=5m, J e = 5 
Fond : 

Cpl e -1750m/s; pl e =l .8; al e =0 .5; hl e =0 . 75m; 

Cp2 e =1970m/s; p2 e =l . 93 ; a2 e -CK4; h2 e =1.75m; 

Cp3 e =1793m/s; p3 e =1.85; a3 e =0.5; h3 e =0.45m; 

Cp4 e =2050m/s; p4 e =2 . ; a4 G =0,38; h4 e =2.25m; 

Cp5 e =3000m/s; p5 e =2 . 3 ; a5 e -0.2; fcs5 e =13 OOm/s ; a5s e =0 . 3 
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Come le comprendra 1 1 homme du metier, au lieu d'utiliser 
une f one ti on de cout telle que K(n) calculee au moyen une 
distance euclidienne, la phase d J analyse ANA pourrait 
tout aussi bien utiliser la recherche d'un maximum de 
5 correlation entre les signaux RfSi n {t, B n ) et Rfi (t,B e ) . 

Comme neanmoins il s'agit, dans 1 f un et 1 ' autre de ces 
deux cas , de reveler des ecarts entre ces signaux et de 
faire evoluer le signal Rf Si n ( t, B n ) dans un sens propre a 
10 reduire ces ecarts, les operations precedemment evoguees 
de soustraction SOUST et de formation de tels ecarts 
doivent etre ici comprises dans un sens fonctionnel large 
incluant notamment une telle recherche de maximum de 
correlation . 



15 
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RE VENDIC ATIONS . 



1. Precede de caracterisation geo technique d'un fond 
5 sub - aqua t i que (B) , tel qu'un fond sous -mar in, recouvert 
d 1 une nappe d 1 eau de hauteur to tale (H) determinee , et 
comprenant une pluralite de strates (Fl, F2) formant 
entre elles des interfaces de separation, presentant 
chacune des parametres physiques propres (hj , pj , C j , aj) 

10 et s 1 etendant a des profondeurs diverses en dessous d'une 
premiere strate (Fl) qui forme une interface de 
separation avec 1 ' eau, ce precede comprenant au mo ins une 
phase d 1 investigation (INVESTIG) , une phase de pre- 
traitement ( PRE_JTRAIT) et une phase d 1 analyse (ANA), la 

15 phase d * investigation ( IWESTIG) incluant elle-meme une 
operation (EXCIT) d 1 excitation acoustique, mise en oeuvre 
en emettant au mo ins vers le fond (B) , depuis une source 
acoustique (1) immergee , une onde acoustique incidente 
(s{t)) de signature f requentielle connue, une operation 

20 d ' acquisition (ACQUI) mise en oeuvre en produisant, au 
moyen d'une flute (2) d'au mo ins quatre recepteurs (21 a 
24) iitimergeS; des signaux de mesure respectifs (Si(t,l)) 
resultant d'une detection d'ondes acoustiques respectives 
ref lechies par le fond (B) , et une operation de releve 

25 (RE LEV) de donnees geometriques (Zs(t), Zr(t), H(t), 
X(t), Y(t)) / la flute (2) etant au mo ins 
approximativement alignee sur la source (1) et distante 
du fond (B) , les recepteurs (21 a 24) etant espaces les 
uns des autres , et 1 1 onde incidente ( s ( t ) ) etant emise 

30 dans chaque fenetre d'une succession de fenetres 
temporelles d* emission disjointes, et presentant une 
variation temporelle de frequence a 1 * interieur de chaque 
fenetre d ' emission, caracterise en ce que la phase 
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d ' investigation (INVE8TIG) inclut une operation de 
deplacement (DEPL) , concomitante avec les operations 
d'excitation (EXCIT) et d 1 acquisition (ACQUI) , et raise en 
ceuvre en deplagant simultanement la source acoustique (1) 

5 et la flute (2) dans la nappe d r eau, en ce que la phase 
de pre-traitement ( PRE_JTRAIT) est mise en ceuvre en 
deduisant , a partir des signaux de mesure si(t,l), de la 
signature de 1 ' onde acoustique incidente (s(t)), et des 
donnees geometriques (Zs { t) , Zr(t), H (t) , X(t), Y(t) } , 

10 des signaux de reflexion rf i ( t, 1) ayant atteint des 
recepteurs differents (21 a 24} a la suite d'une m§rne 
emission d 1 onde incidente, et en regroupant , sous forme 
de series de signaux de reflexion recales (Rf i { t , B) ) , des 
signaux de reflexion rfi(t,l) ayant atteint, sous des 

15 angles de rasance (6±) differents, des recepteurs 
differents (21 a 24) en provenance d'une meme zone du 
fond a la suite de dif f erentes emissions de 1 * onde 
incidente, et en ce que la phase d 1 analyse (ANA) est mise 
en oeuvre en deduisant les parametres physiques {hj , pj , 

20 Cj , cxj ) de chaque strate (Fl a F4) du fond, a partir des 
series de signaux de reflexion recales (Rf i { t , B) ) , par 
exploitation de la dependance angulaire et f requentielle 
du coefficient de reflexion de 1 1 onde incidente sur 
chaque interface, chaque signal de reflexion recale 

25 <Rfi<t,B)) etant interprets comma une somrae d'arrivees 
retardees et attenuees de rayons propres correspondant a 
la reflexion de 1 1 onde incidente par les dif f erentes 
strates de la zone du fond a laquelle correspond ce 
signal de reflexion recale (Rfi ( t , B) ) . 

30 

2. Procede de caracterisation geo technique suivant la 
revendication 1, caracterise en ce que chaque onde 
ref lechie atteignant un recepteur est issue de la 
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reflexion d'une onde incidente ayant atteint le fond sous 
un angle de rasance {9±) au plus egal a 7 0 degres et de 
preference au plus egal a 65 degres. 

5 3 . Procede de caracterisation geotechnique suivant 
l'une quelconque des revendications precedentes, 
caracterise en ce que chaque onde ref lechie atteignant un 
recepteur est issue de la reflexion d'une onde incidente 
ayant atteint le fond sous un angle de rasance (Q x ) au 
10 moins egal a 10 degres et de preference au mo ins egal a 
2 5 degres. 

4 . Procede de caracterisation geotechnique suivant 
l'une quelconque des revendications precedentes, 

15 caracterise en ce que chaque recepteur (21 a 24) est 
separe d 1 un recepteur voisin (22 a 24) par une distance 
(d) • correspondant a une difference d' angle de rasance au 
moins egale a 2 degres, et de preference au moins egale a 
5 degres, des signaux de reflexion respectivement regus 

20 par ce recepteur et par le recepteur voisin. 

5 . Procede de caracterisation geotechnique suivant 
l'une quelconque des revendications precedentes, 
caracterise en ce que 1 1 onde incidente presente, dans 

25 chaque fenetre d' emission, au moins une composante de 
frequence relativement basse, dont la frequence est au 
moins egale a 10 0 Hz. 

6 . Procede de caracterisation geotechnique suivant 
30 l'une quelconque des revendications precedentes, 

caracterise en ce que 1 1 onde incidente presente, dans 
chaque fenetre d 1 emission, au moins une composante de 



10 
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frequence relativement haute, dont la frequence est au 
plus egale a 8 kHz. 

7. Precede de caracterisation geotechnique suivant 
1'une quelconque des revendications precedentes , 
caracterise en ce que 1 - onde incidente presente, dans 
chaque fenetre d' emission, une frequence contingent 
variable entre une composante de frequence relativement 
basse et une composante de frequence relativement haute. 




8. Precede de caracterisation geotechnique suivant 
1'une quelconque des revendications precedentes, 
caracterise en ce que les fenetres temporelles d ' emission 
ont chacune une duree au moins egale a 0.1 seconde et de 

15 preference au moins egale a 0.5 seconde. 

9. Precede de caracterisation geotechnique suivant 
1'une quelconque des revendications precedentes, 
caracterise en ce que les fenetres temporelles d- emission 

20 sont repetees periodiquement avec une periode au moins 
egale a 0.5 seconde et de preference au moins egale a 1 
seconde . 

10. Procede de caracterisation geotechnique suivant 
25 1'une quelconque des revendications precedentes, 

caracterise en ce que les parametres physiques sont au 
moins en partie choisis dans 1- ensemble de parametres 
comprenant 1 ' epaisseur (h) , la densite (p) , l a vitesse du 
son {C p , C.), l'attenuation du son (« p> « s ) , la rugosite, 
30 et un gradient de celerite du son. 



11. Procede de caracterisation geotechnique suivant 
1'une quelconque des revendications precedentes, 
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caracterise en ce qu'il comprend, en amont de la phase 
d ' analyse (ANA) , une operation (CAPT) de determination de 
la prof ondeur d' immersion (Zs) de la source acoustique 
et/ou de la prof ondeur d' immersion (Zr) d'un point de la 
5 flute et/ou de la hauteur d'eau H. 

12 . Precede de caracterisation geotechnique suivant la 
r e vendi c a t i on 10 , caracterise en ce que la phase 
d* analyse (A3STA) est mise en (Euvre en parcourant, de fagon 

10 iterative, une boucle incluant une operation de 
modelisation (MODEL_FOND) au cours de laquelle sont 
calcules, a partir de parametres physiques fictifs (h nj/ 
Pnj, C nj/ Q'nj ) virtue llement attribues a la premiere strate 
(Fl) au moins de la pluralite de strates, et pour les 

15 differents recepteurs (21 a 24), des signaux de reflexion 
recales fictifs respectifs (Rfs i n (t / B n ) ) issus de signaux 
de mesure virtue llement produits par ces recepteurs (21 a 
24) en presence de ces parametres fictifs (h nj , p nj , C n j , 
^nj ) / une operation de soustraction (SOUST) au cours de 

20 laquelle sont formes des ecarts (A ni ) constates entre les 
signaux de reflexion recales respectifs (Rfi(t,B)) issus 
des signaux de mesure (si(t f 1)) reellement produits par 
ces recepteurs (21 a 24) et les signaux de reflexion 
recales fictifs correspondants (Rfs^Ct, B n ) ) , et une 

25 operation d 1 optimisation et de retroaction (OPTIJRETR) a 
1 1 issue de laquelle les parametres physiques fictifs sont 
revises en fonction des ecarts constates. 
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